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ABSTRACT 

Until 2004, anomalous X-ray pulsars (AXPs) were known as strong emitters of soft X-rays only (< 
10 keV). The discovery of hard X-ray component from AXPs provided important insight about their 
emission properties while it posed a serious challenge to explain its origin. The physical mechanism 
of the hard emission component has still not been fully resolved. We investigate the high-energy 
gamma-ray properties of the brightest AXP, 4U 0142+61 using data collected with the Large Area 
Telescope on board Fermi Gamma-ray Space Telescope to establish the spectral behavior of the source 
on a very broad energy span and search for pulsed emission. Here, we present our results of detailed 
search for the persistent and pulsed high-energy gamma-ray emission from 4U 0142+61 which result 
in no significant detection. However, we obtain upper limits to the persistent high-energy gamma-ray 
emission flux which helps us to constrain existing physical models. 
Subject headings: pulsars: individual (AXP 4U 0142+61 ) — gamma-rays: stars 



1. INTRODUCTION 

Anomalous X-ray pulsars (AXPs) have been in- 
triguing sources since the ir discovery in early 1980s 
(jFahlman fc Gregory! Il98l . They are bright X-ray 
sources with X-ray luminosities (below 10 keV) in the 
range of 10 33 -10 36 erg s _1 . They spin rather slowly and 
their spin periods clustered in a narrow range of 2-12 
s. The lack of any evidence for binary nature (e.g., 
Doppler modulation in their long term pulse periods) 
eliminates the possibility that they are accreting mat- 
ter from a donor. Their spin down rates are relatively 
large, ranging between 10 -10 and 10~ 12 s s _1 . Their 
rotational energy loss is insufficient by orders of magni- 
tudes to provide the observed X-ray luminosities. AXPs 
are commonly regarded as young isolated neutron stars 
that are powered b y their extremely strong m agnetic 
fields, B> 10 14 G (jPuncan fc Thompson! IT992T) . Such 
strong magnetic fields can efficiently slow these young 
systems down via magnetic breaking and provide energy 
for the em itted X-rays via diffusion o f evolving mag- 
netic field (jThompson fc Duncan! 1 1 9961). A detaile d re- 
view on AXPs can bejbund in iMereghettH (|2008) and 
iWoods fc Thompson! (|2006t ). 

A major observational development in AXP 
studies was the discovery of hard X-ray emis- 
sion from AXPs IE 1841-045 (IKuiper et al.l 
12001 . 4U 

1RXS J170849. 0-400910 



0142+61 (Iden Hartog et al.l 120041) an d 
(jRevnivtsev et al.l 12004!) . 
IE 1841-045, located at the center of supernova remnant 
(SNR) Kes 73, is the first AXP from which non-thermal 
pulsed hard X-ray /soft ga mma-ray emission was dis- 
covered (jKuiper et al.l 12004). The pulsed nature of this 
emission eliminates the possibility of its SNR origin. 
Spectral studies of this source using INTEGRAL obser- 
vations in the 20—300 keV band revealed a power— law 
shape with an i ndex of 1.39 ± 0.05 but no evidence of 
spectral break (jKuiper et al.l [2006) . The hard X-ray 



1 Sabanci University, FENS, Orhanh- Tuzla, Istanbul 34956 
Turkey 



emission from 1RXS J170849. 0-400910 was discovered 
in the Galactic Plane Survey observations with INTE- 
GRAL dRevnivtsev et a l. 2004) and shown to be puls ed 
as well (jKuiper et al.l (|2006l) . Iden Hartog et al.l (|2008f) V 
Finally, hard emission component from 4U 0142+61 was 
discovered during INTEGRAL observati ons of the 
Cassiopeia region (jden Hartog et al.l 120041) . Detailed 
studies with 2.37 Ms INTEGRAL observations showed 
a power law spectrum up to about 230 keV with an 
index of 0.93 ± 0.06 (jden Hartog et al.ll2008l ). Based 
on the logparabolic function fit to the INTEGRAL SPI 
and ISGRI observations, they estimate a peak energy 
of the spectral energy distribution to be ^228 keV 
(jden Hartog et alj 12008). and 20-150 keV flux as (8.97 
± 0.86) x 10 -11 erg cm" 2 s . It is important to note 
that the energy emitted above 15 keV is comparable or 
larger than that in the soft X-ray band (that is, below 
10 keV, see e.g.. Iden Hartog et all (|2008f )). 

The physical nature of the hard X-ray / soft 
gamma-ray emission is s till n ot well understood. 
IBeloborodov fc Thompson! (|2007f) proposed that the 
hard X-ray emission could originate from a plasma 
corona around the magnetar. They suggest that such a 
corona around magnetars can be formed via starquakes 
that could she ar the neutron star crust and its external 
magnetic field. IHevl fc Hernquistl (|2005al Tbh proposed the 
fast-mode break down model in which an optically thick 
fireball produced by the magnetohydrodynamics waves 
created near the surface of the neutron star. Further 
they suggest that if fast modes are not strong enough 
to yield an optically thick fireball, the produced non- 
thermal emission would be sufficient to explain the ob- 
served high-energy emission from soft gamma repeaters 
and AXPs. A nother attempt to explain the hard X-ray 
emission is bv iBaring fc Harding! ((200l . They suggest 
that the upscattered surface thermal X-ray photons as 
the source of observed non-thermal hard X-ray emission 
from AXPs. 

In order to understand the nature of the hard emis- 
sion component of AXPs, it is crucial to establish their 
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spectral shapes on a wide range in the energy domain. 
In particular, it is important to determine where their 
spectral energy distribution peak. Thanks to the Large 
Area Telescope (LAT) on board Fermi Gamma-ray Space 
Telescope (Fermi) , we are now able to investigate high- 
energy behavior of these sources with an unprecedented 
data quality. 

In this paper, we performed detailed search for persis- 
tent and pulsed high-energy gamma-ray emission from 
4U 0142+61 using Fermi/LAT observations. We also 
employed contemporaneous Rossi X-ray Timing Explorer 
(RXTE) observations to obtain the spin ephemeris of the 
source. In Section 2, we describe the observations used 
and details of data analysis. We present our results in 
Section 3 and discuss their implications in Section 4. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Fermi/LAT 

LAT is one of the two instruments on board Fermi 
operating in the energy band of 20 MeV— 300 GeV. 
It is a pair conversion telescope with a high-resolution 
silicon tracker, calorimeter, anticoincidence detector, 
programmable trigge r and data acquisition system 
(|Atwood et~ani2009h . Since 2008 August 4, the LAT 
has been operating as an all sky monitor in high-energy 
gamma-rays, covering the full sky in approximately every 
3 hr. 

We accumulated the LAT data within a 15° radiusQ 
centered at 4U 0142+61, collected from 2008 August 4 
to 2010 April 29 with an exposure time of ~31.7 Ms. 
We also performed spectral analysis using a 2° radius 
region around the source in order to completely avoid 
contamination from the nearby bright sources . We per- 
formed our unbinned likelihood analysis using Science- 
Tools v9rl5p2 with P6_V3_DIFFUSE set as the instru- 
mental response. In the event selection process, we set 
the maximum zenith angle to 105° in order to eliminate 
background gamma-rays due to the Earth limb. All time 
intervals where the zenith cut intersects the region were 
excluded. In Figure 12.11 we present the count map im- 
age of the 40° radius region centered at 4U 0142+61 . 
The diffuse gamma-ray emission from the Milky Way was 
modeled with the latest model, gll_iem_v02.fit. We also 
used isotropic_iem_v02.txt to account for the extragalac- 
tic isotropic diffuse emission and residual instrumental 
background. The spectral fits and flux calculations were 
done with the python version of gtlike, pyLikelihood. 

For timing analysis, photon arrival times of all events 
within the 2° extraction region around 4U 0142+61 were 
converted to that at the solar system barycenter using 
gtbary of ScienceTools. 

2.2. RXTE 

To search for pulsed high-energy gamma-ray emission 
from 4U 0142+61, we obtained the precise spin ephemeris 
of the source using contemporaneous RXTE observa- 
tions. On board RXTE, there are three scientific pay- 
loads: the Proportional Counter Array (PCA) that is 
sensitive to photon energies between 2—60 keV, the High 
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Energy X-ray Timing Experiment, sensitive to photons 
in the 15-250 keV photons and the All Sky Monitor. We 
have employed only the PCA observations to achieve our 
goal. 

4U 0142+61 has been monitored with the RXTE pe- 
riodically for the last ^8 years with pointings almost 
uniformly spaced, usually by about two weeks. We have 
selected 53 RXTE observations that were performed be- 
tween 2008 August 4 and 2010 April 30 (under the Pro- 
gram IDs: P93019, P94019 and P05019) which cover 
the investigated LAT observing span. Individual RXTE 
pointings are typically between 3 and 4 ks long and the 
total exposure time of all selected observations is about 
196 ks. For each observation, we extracted events in the 
2—10 keV range collected with the PCA and converted 
their arrival times to the solar system barycenter using 
the faxbary tool of HEASoft 6.8.' 

3. RESULTS 
3.1. Search for Persistent Emission 

As evident from Figure [2~T1 4U 0142+61 is clearly not 
detected in the LAT energy passband. A point source 
search using the filtered event list with the gtfindsrc tool 
of ScienceTools results in a potential source whose coor- 
dinates are inconsistent with that of 4U 0142+61, there- 
fore yields no detection. 

In order to obtain very-high-energy gamma-ray flux 
upper limits of 4U 0142+61, we fitted the data from the 
15° radius region. We added all br ight cataloged sources 
and recently discovered blazar ( Vandcnbr oucke et al.l 
2010) within this region of interest into the model as well 
as the galactic diffuse and extragalactic diffuse emission 
leaving their model parameters free. The fit yields a test 
statistics (TS) value of ~0.23 which implies a detection 
significance less than la. The 3<r flux upper limits with 
a power law index 2.5 are 2.32 x 10~ 6 MeV cm~ 2 s _1 in 
the 0.2-1.0 GeV band and 1.28 x 10~ 6 MeV cm" 2 s" 1 
in the 1.0-10.0 GeV band. Note here that the spectral 
parameters and fluxes of the cataloged sources obtained 
in the latter fit are consistent with the catalog valued, 
showing that our analysis is robust. 

We performed similar spectral modelling for the data 
of the 2° radius region with a power law model of index 3 
as well as the galactic diffuse and extragalactic isotropic 
diffuse emission models. The resulting TS value is ^3 
which implies a detection significance less than 2a. We 
chose the 0.2-1.0 and 1.0-10.0 GeV energy bands for flux 
calculations and find 3<r upper limits to the source flux 
in these energy bands as 5.72 x 10~ 6 MeV cm~ 2 s" 1 and 
1.29 x 10~ 6 MeV cm~ 2 s -1 , respectively. In Figure |H 
we present the high-energy gamma-ray flux upper limits 
of 4U 0142+61 in the vF v repre sentation along with its 
low energy gamma-ray behavior (den Har tog et al.ll2008l 
the data obtained from). We discuss their implications 
in §4. 

3.2. Search for Pulsed Emission 

We employed a Fourier based epoch folding tech- 
nique to obtain the spin ephemeris of 4U 0142+61 us- 
ing RXTE /PCA observations covering the time span of 
the LAT exposure of the source. We first generated the 
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Fig. 1. — Smoothed Fermi-LAT count map around AXP 4U 0142+61 in the 0.2—100 GeV energy range. The large circle and small circle 
show the 15° and 2° radius extraction regions, respectively. The plus sign indicates the position of AXP 4U 0142+61. 



pulse profile of the source using three consecutive PCA 
observations around the epoch (MJD 54713.5). Then, 
we grouped observations in order for them to be spaced 
at least 0.2 days apart from each other, and we cross 
correlated the pulse profile of each group of pointings 
with the template profile to determine the phase shift of 
each pointing with respect to the template. Finally, we 
fit the phase shifts with a polynomial to obtain the spin 
ephemeris. In Figure [3~2| we present the phase shift and 
the best fitting model, that is a third-order polynomial 
(x 2 /degrees of freedom = 59.6/42). We tabulate the best 
fit spin ephemeris parameters of 4U 0142+61 in Table [TJ 



TABLE 1 
Spin ephemeris of 4U 0142+61 



Parameter 


Value 


Range (MJD) 


54682.6 - 55315.1 


Epoch (MJD) 


54713.5 


v (Hz) 


0.1150900026(9) 


V (10~ 14 Hz s" 1 ) 


-2.745(8) 


V (10~ 23 Hz s" 2 ) 


3.6(3) 



To search for pulsed high-energy gamma-ray emission 
from 4U 0142+61, we generated the LAT pulse profiles 
of the source in the 0.2-1.0 GeV and 1.0-10.0 GeV energy 
ranges using the precise PCA spin ephemeris obtained. 
We find that both LAT profiles are consistent with ran- 
dom fluctuations with respect to its mean. We calculate 
a 3<7 upper limit to the rms pulsed fraction of 1.5% in the 
0.2-1.0 GeV band and 2.3% in the 1.0-10.0 GeV band. 
We also investigated the lower energy part of the LAT 
passband (30-200 MeV), which also resulted with no ev- 
idence of pulsed emission; the 3<r rms pulsed fraction 



upper limit is 1.6%. 

4. DISCUSSION 

We searched for both persistent and pulsed high-energy 
emission from the AXP 4U 0142+61 using Fermi /LAT 
data. We find no significant detection in either of the two 
objectives. Nevertheless, we obtained 3<r upper limits for 
the high energy persistent emission in the 0.2-1.0 GeV 
and 1.0-10.0 GeV ranges of 5.72 x 10" 6 MeV cm' 2 s' 1 
and 1.29 x 10 -6 MeV cm -2 s _1 , respectively. As for 
the pulsed emission, a 3u upper limit to the RMS pulsed 
amplitude in the 0.2-1.0 GeV range is 1.5% and in the 
1.0-10.0 GeV range is 2.3%. The search in the lower 
energy LAT passband (30-200 MeV) also did not yield 
a pulsed emission. The 3<r rms pulsed amplitude upper 
limit is 1.6%. Our LAT upper limits are much lower than 
high pulsed fractio n (up to 100%) seen in hard X-rays 
with INTEGRAL (jden Hartog et al.ll2008h . 

In order to establish the spectral shape of 
4U 0142+61 on a wide energy range, we constructed 
the v¥ v spectrum of the source in the 15 keV— 10 GeV 
ra nge by adopting the ha rd X-ray spectrum presented 
in Iden Hartog et al.1 (l2008|) and placing the upper limits 
calculated in this work. " Iden Hartog et al.1 (|2008| ) fitted 
the INTEGRAL /ISGRI data with a simple power 
law model of index 0.93 ± 0.06. We present this fit 
with dashed lines in Figure |4j We place an upper 
limit curve, which is the line connecting the two LAT 
upper limit measurements (as shown with solid line in 
Figured]), is also a power law with an index of —0.76. 
These two curves intersect with each other at ~1.1 
MeV which is an upper limit to the spectral break 
energy. Note that the spectral brea k upper limit is 
consistent with Iden Hartog et al.l (|2008|) measurement of 
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Fig. 2. — (Top) spin phase shifts of RXTE/PCA observations of 4U 0142+61 with respect to the epoch. The solid line is the best fitting 
model, that is a third-order polynomial. (Bottom) residuals of the fit. 



2791^ keV obtained by fitting a logparabolic function 
to the combined XMM-Newton, INTEGRAL /ISGRl, 
INTEGRAL/SPl and CGRO /COMPTEL data. 

Our estimated upper limit to the spectral break en- 
er gy is in accordance with the coro nal emission model 
bv lBeloborodov fe ThompsonI ([2007). According to their 
model, photons with energies in excess of >1 MeV 
would be trapped in the ultrastrong magnetic fields 
(B > 10 14 G). In such a case photons would cither 
split into two photons or they would create an electron- 
positron pair, therefore, suppressing the emission from 
the inner corona above ~1 MeV. Our estimate of the 
spectral break energy upper limit also agrees with the 
predictions of the quantu m electrodynamics model by 
iHevl fc HernquisH (2005a. b) as they expect the lower 
limit to the break energy to be around 1 MeV. On the 
other hand, they claim that if a source has a signifi- 
cant excess emissio n in optical waveleng ths, as in the 
case of 4U 0142+61 (jHulleman et al.ll2000l ). its v¥ v spec- 



trum should continuously increase in the 10-200 MeV 
range. Note that optical emission originates from the 
neutr on star itself (|Kern fe Martini [2002; Dh illon et al. | 
2005[) and from a disk around 4U 0142+61 (Wang et all 
20061 lErtan fe Qahskanl [20061 ). However, the origin of 



the excess optical emission in 4U 0142+61 is not clear. 
If the excess optical emission originates from the com- 
pact object, our results place constraint, although not 
stringent, on the quantum electrodynamics model due 
to the lack of increase in the 10-200 MeV range in the 
v¥ v spectrum. On the other hand, if the disk provides 
a significant contribution to the observed emission, then 
the optical radiation from the neutron star itself would 
not be excessive and the quantum electrodynamics model 
would still remain feasible. 



We thank M. Ali Alpar and Yuki Kaneko for valu- 
able discussions. S.ij.M. and E.G. acknowledge EU FP6 
Transfer of Knowledge Project Astrophysics of Neutron 
Stars (MTKD-CT-2006-042722). 



REFERENCES 



Atwood, W. B., et. al. 2009, ApJ, 697, 1071 
Baring M. G., & Harding, A. K. 2007, Ap&SS, 308, 109 
Beloborodov, A. M. & Thompson, C. 2007, ApJ, 657, 967 
den Hartog, R. P., Kuiper, L., Hermsen, W., & Vink, J. 2004, 
ATel, 293, 1 

den Hartog, R. P., Kuiper, L., Hermsen, W., Kaspi, V. M., Dib, 
R., Knodlseder, J., & Gavriil, F. P. 2008, A & A, 489, 245 

Dhillion, V. S., Marsh, T. R., Hulleman, F., van Kerkwijk, M. H. 
Shearer, A., Littlefair, S. P., Gavriil, F. P., & Kaspi, V. M. 
2005, MNRAS, 363, 609 

Duncan, R. C, & Thompson C. 1992, ApJ, 392, L9 

Ertan, U., & Qahskan, S. 2006, ApJ, 649, L87 



Fahlman, 
Heyl, J. S 
Heyl, J. S 
Hulleman. 

408, 68! 
Kern, B., 
Kuiper, L 
Kuiper, L 



G. G., & Gregory, P. C, 1981, Nature, 293, 202 

., & Hernquist, L. 2005a, ApJ, 618, 463 

., & Hernquist, L. 2005b, MNRAS, 362, 777 

, F., van Kerkwijk, M.H., Kulkarni, S.R., 2000, Nature, 



l Martin, C. 2002, Nature, 417, 527 
Hermsen, W., Mendez, M. 2004, ApJ, 613, 1173 
Hermsen, W., den Hartog, P. R. 2006, ApJ, 645, 556 
Mereghetti, S., 2008, A&AR, 15,225 
Revnivtsev, M. G., et al. 2004, Astron. Lett., 30, 382 
Thompson, C, & Duncan, R. C. 1996, ApJ, 473, 322 
Vandenbroucke, J. et al. 2010, ApJ, 718, L166 



10" 



10~ 2 10" 1 10° 10 1 10 2 10 3 10 4 

Energy [MeV] 

Fig. 3.— Wide band vF„ spectrum of 4U 0142+61: INTEGRAL/lSGRl (20-300 keV) in black (stars), INT EGRAL /SPY (20-1000 keV ) 
in red (open squares) and CCRO/COMPTEL (0.75-30 MeV) 2<r upper limits in black (data constructed from lden Hartog et all {2008)). 
Fermi/LAT upper limits (in the 0.2-1.0 GeV and 1.0-10.0 GeV) obtained using the 2° extraction region are shown in b lue diamonds and 
that of 15° extraction region in red triangles. Dashed line is the best fit power law model to the ISGRI data points (den Hartog ct al. 
2008). Solid line shows the power law upper limit trend of the 2° Fermi/LAT region. 
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